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Abstract. A classical approximation to time dependent quantum mechanical 
scattering in the M0ller formalism is presented. Numerically, our approach is 
similar to a standard Classical— Trajectory— Monte— Carlo calculation. Conceptu- 
ally, however, our formulation allows one to release the restriction to stationary 
initial distributions. This is achieved by a classical forward-backward propagation 
technique. As a first application and for comparison with experiment we present 
fully differential cross sections for electron impact ionization of atomic hydrogen 
in the Erhardt geometry. 



PACS numbers: 34.80D, 03.65. Sq, 34.10+x 

Classical models and approximations are frequently used for atomic and molecular 
problems, despite their inherent quantum nature. One of the reasons is that our 
intuitive understanding is mainly based on classical terms and pictures by which 
we are surrounded in every day life. Another reason for not doing (fully) quantum 
mechanical calculations is the complexity of a problem: fully differential cross sections 
in higher dimensional atomic systems, e.g., often require a numerical effort still beyond 
present computing power. Only recently the quantum mechanical Coulomb three- 
body scattering problem was solved numerically (l], ^| . 

On the other hand, a remarkably successful classical approach to collisional atomic 
problems has been developed over the years, the so called Classical-Trajectory-Monte- 
Carlo method (CTMC). It was introduced as a purely classical model based on a 
"planetary atom" with a major axis of two meters (!) ||. This model has produced 
reasonable results for total or energy differential ionization cross sections on atomic 
hydrogen || and for other few-body Coulomb collision processes Q. 

Attempts to reduce the limitation of this classical model aimed at changing 
the description of the initial state for the hydrogen atom from a microcanonical 
distribution to one, that is closer to the quantum density |7[ |). Another idea was 
to introduce additional ad hoc stabilisation potentials in order to be able to treat 
multi-electron targets ||. However, all attempts took as a starting point not the 
quantum problem but the previously formulated classical model. Hence, the proposed 
amendments were accompanied by inconsistencies or the need of "fit parameters" 
determined from cross sections. In the end, one must say that it is up to now not 
possible to describe higher differential cross sections or targets with more than one 
active electron consistently in a classical collision framework. 
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To achieve progress in this situation we decided to go one step backwards and 
start with a time dependent quantum mechanical scattering formalism. By following 
the approximations which lead to the classical description, i.e. the CTMC method, we 
can identify the source and nature of deviations between the classical and the quantum 
result which serve as a guide to improve the classical description to a quasi-classical 
approximation. 

We divide the problem into three logically separate steps: (1) preparation of 
the initial state before the collision, (2) propagation in time, and (3), extraction of 
the cross section. For a consistent quasiclassical picture each of these steps has to be 
approximated in the same way. To keep the derivation transparent, we will concentrate 
on electron impact ionization of one active target electron in the following. 

Step (1): The initial wave function for the collision problem is translated first 
into a quantum phase space distribution by the Wigner transformation (hJ. The 
resulting Wigner distribution is reduced to a classical distribution w(p, q) which can be 
propagated classically in phase space by taking the usual limit T% — > 0. The difference 
to an a priori classical approach is the use and interpretation of negative parts of the 
distribution: Viewed as the fi, — > limit of a quantum problem they contribute to 
the observables in the same way as the positive parts since they do not need to be 
interpreted as weights for real paths of classical particles. Yet, there arise additional 
problems when using this type of general initial phase space distributions in the usual 
classical framework: Most of them are not stationary under classical propagation, their 
Poisson bracket with the Hamilton function does not vanish, {H, w} ^ 0. Hence, 
the initial target distribution will look very different at the time the projectile has 
approached and the collision actually happens. 

Step (2): The formulation of the propagation is crucial since it must resolve 
the problem of the non-stationary classical initial distribution, as described above. 
Traditionally, the time dependent scattering is described by calculating the transition 
amplitude between initial and final state through the S-matrix, which is in turn related 
to the t-matrix describing directly the cross section, see, e.g., |0. In a simplified 
version where the asymptotic initial and final states are eigenstates of the asymptotic 
Hamiltonians and Hq one normally writes for the transition amplitude 

S fi = lim (f\U(t)\i) (1) 

where U(t) = exp[— iHt] denotes propagation with the full Hamiltonian. By a Wigner 
transform the quantum time evolution operator U(t) can be directly transformed with 
the help of the quantum Liouville operator C qi which reduces to the classical Liouvillc 
operator C c in the limit h — ► The latter describes the evolution of a phase 

space distribution u>(r, p, t) according to the Poisson bracket 

d t w = {H, w} = —i£ c w (2) 

in analogy to the quantum evolution of the density matrix p generated by the 
commutator, 

d t p=-i[H,p}. (3) 

Hence, we could directly use the translation of ([!]) to classical mechanics via the 
Liouville operator. In connection with the microcanonical initial state distribution this 
is indeed equivalent to the CTMC formulation ]l3| . However, using non-stationary 
initial state distributions is inconsistent with the reduced quantum description of (m) 
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which relies on the fact that the asymptotic states are eigenstates of Uq and therefore 
stationary. Instead we have to go back to the full scattering formulation 

S fi = (f\SllSl+\i), (4) 

where 

T = lim U\t)U (t) (5) 

f— >±oo 

are the M0ller operators. The meaning of fi+, e.g., is to propagate backwards with 
Uo(t) using the asymptotic Hamiltonian without the projectile-target interaction 
and then forward again under the full Hamiltonian with U (t) . Again, with the help of 
the Liouville operator we can translate the M0ller operators to their classical analogue, 
thereby obtaining a prescription how to propagate a non-stationary initial phase space 
distributions ^(7), where 7 = (pi,qi,p 2 ,q 2 ) is a point in the 12-dimensional phase 
space: 

w f = lim lim e- iC - t e iC " ( - t - t ">e iC - t 'wi=K.Wi. (6) 

t— >+oo t'^ — 00 

The difference to ([!]) are the explicit propagations under C{. and C l c in the 
initial and final channel (which need not be the same). The meaning of (^) becomes 
very transparent if we insert a discretized distribution, which is used in the actual 
calculations, 1(^(7) = J2n w n^ 12 {l ~ l n )- The final distribution reads 

Wf (7) = K-Wi = WnSl2 (T - 7n) (?) 

n 

where each phase space point 7^ emerges from 7^ through solving successively 
Hamilton's equations, first with Hq, then with H, and eventually with . With 
this propagation scheme a non-stationary initial distribution will spread when being 
propagated backwards with the asymptotic C l c . However, it will be refocused under the 
following forward propagation with C c . Hence, when the actual collision happens for 
t « the original target distribution is restored, slightly polarized by the approaching 
projectile. 

Hence, there is no more need for the initial distribution to be classically stationary. 
We are able to use any phase space distribution as a target in our quasi classical 
approach. This also includes unstable multi-electron targets, e.g., classical helium. 

Step (3): Before we come to the actual evaluation we have to formulate the 
cross section such that it can make full use of the non-stationary initial phase space 
distribution Wi(pi,qi), where "I" refers to the target electron. Without modification 
the total energy E of the final state forces by energy conservation for each classical 
trajectory only those parts of the initial phase space distribution to contribute to the 
cross section which have the same energy E. However, this would bring us essentially 
back to the microcanonical description. In order to make the entire non-stationary 
initial state distribution "visible" to the collision process, we use the energy transfer 
Ei = E± —Ei to the target electron rather than its energy E[ f) itself as a differential 
measure. Of course, as long as the initial state is on the energy shell with well defined 

(i) (i) 

energy E = EY> + E£' the new definition coincides with the usual expression for the 
cross section, 

(8) 



dtt 1 dtt 2 dE 1 



dVLidVL 2 dE 1 



where dVti are the differentials for the solid angles of the two electrons, respectively. 
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To extract this cross section we have to evaluate the phase space integral 



of 



dtt 1 dtt 2 dE 1 



Jdx 2 dy 2 dp, dq, f[ 6{nP - ni)S(E[ f) - E 1 )w i , (9) 

i—l 

where the integration is over the initial state variables, namely the impact parameter 
area dx 2 dy 2 and the phase space of the (bound) target electron dpidqi, with initial 
distribution Wi(jpi, q\, x 2 , y 2 ). The propagated angles £l{ of the electrons as well as the 
energy transfer El have to coincide with the desired values f2i and E\ to contribute to 
the cross section (|)|) which is a generalization of the one derived in e.g., where the 
initial bound state was assumed to live on a torus, i.e., Wi{p\,q{) = 5(I(pi,q\) — Iq) 
with a well defined multidimensional action Iq. 

Finally, we have to respect the Pauli principle for the two identical electrons. 
Formally, this can be done easily in the Wigner transform for the two ionized 
electrons in the final state. In the limit Ti — one is left with the usual classical 
symmetrization, i.e., an interchange of indices. To keep the notation simple we have 
omitted symmetrization in the outlined derivation, however, it is included in the actual 
computation which is carried out by applying standard CTMC techniques to evaluate 
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Figure 1. Fully differential cro^ section in Ehrhardt geometry: Ei„ = 
250eV, Eb= 5eV and 9 a = 8°: comparison of measurements of Erhardt et al 
to quasiclassical calculations with (a), the microcanonical distribution (standard 
CTMC), and (b) with the Wigner distribution. Error bars indicate the statistical 
error of the Monte Carlo calculations. The theoretical data has been scaled by 
about 20% to reproduce the correct total ionization cross section at 250 eV. For 
the negative parts in the Wigner cross section see text. 



As a first application we discuss fully differential ionization cross sections 
for atomic hydrogen comparing three data sets: a calculation with the 
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Figure 2. Same as in Figure 1, but i3r£L n = 54.4eV, E b = 5eV and a = 23° and 
measurements (circles) by Roder et al lid] . 



standard microcanonical distribution (CTMC), one with the non-stationary Wigner 
distribution in our quasi classical framework and experimental data at impact energies 
of 250eV (fig. % and 54.4eV (fig. |), respectively. 

For each impact energy about 10 8 trajectories have been calculated. The cross 
section at 250 eV for the Wigner distribution still exhibits negative parts, indicating 
that this cross section is not yet fully converged. This is not surprising if one takes 
into account that the fraction of phase space of the final state is so small with the 
chosen bin sizes for energies and angles, that only between 100 and 300 events finally 
contribute to the shown cross sections. However, a considerable advantage of the 
present method is that a sampling of 10 8 trajectories contains the complete scattering 
information, not just one specific differential cross section. 

The figures show, that the microcanonical distribution, i.e. standard CTMC, is 
not able to reproduce the binary peak flffl , whereas with the Wigner distribution it is 
reproduced fairly well for 250eV and rather well for 54.4eV impact energy. Keeping 
in mind that in contrast to the microcanonical distribution the Wigner distribution 
has the correct probability densities in momentum and configuration space, one can 
conclude, that at least for energies between 50 and 250 eV the differential cross sections 
"image" the initial phase space distribution. The present approach is still a classical 
approximation and cannot reproduce quantum effects. Therefore, features in the cross 
section, for which coherence is crucial are represented purely or not at all. 

To summarize, we have shown, that a consequent classical approximation to a 
quantum system can give much better results compared to those from an a priori 
classical model, though both approaches are realized numerically in almost the same 
way. However, the main difference is conceptual: in the usual classical limit each 
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individual trajectory represents that of an electron obeying the classical equations of 
motion, whereas in our classical approximation only the entire phase space distribution 
is meaningful and individual trajectories are only discretized points of the distribution 
evolving in time. Hence, there is no problem to deal with "negative probabilities" in 
the initial distribution, since we regard them not as probabilities but only as weights 
of the integration which, of course, may be negative. 

However, the use of non-stationary distributions like the Wigner distribution as 
an initial state implies additional difficulties for a scattering description which we 
have overcome by using a forward-backward propagation scheme akin to the quantum 
M0ller formalism and by a reformulation of the energy differential cross section in 
terms of the energy transfer during the ionization process. With these modifications 
all the tools of the standard CTMC technique can be applied straightforwardly. 

Moreover, our approach can be in principle generalized to multi-electron targets 
since we generate our initial phase space distribution from a quantum wave function 
and we know how to deal with non-stationary initial distributions. 

Financial support by the Deutsche Forschungsgemeinschaft within the SFB 276 
at the University Freiburg is gratefully acknowledged. 
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